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ABSTRACT 



Earth-mass dark matter microhalos with size of ~ 100 AUs are the first struc- 
tures formed in the universe, if the dark matter of the Universe are made of 
neutralino. Here, we report the results of ultra-high-resolution simulations of 
the formation and evolution of these microhalos. We found that microhalos have 
the central density cusps of the form p oc r~^'^, much steeper than the cusps 
of larger dark halos. The central regions of these microhalos survive the en- 
counters with stars except in very inner region of the galaxy down to the radius 
of a few hundreds pes from the galactic center. The annihilation signals from 
nearest microhalos are observed as gamma-ray point-sources (radius less than 
1'), with unusually large proper motions of ~ 0.2 degree per year. Their surface 
brightnesses are ~10% of that of the galactic center. Their S/N ratios might 
be better if they are far from the galactic plane. Luminosities of subhalos are 
determined only by their mass, a nd they are mo r e than one order of magnitude 
luminous than the estimation by ISpringel et al.l ( 120081 ): A boost factor can be 
larger than 1000. Perturbations to the millisecond pulsars by gravitational at- 
tractions of nearby earth-mass microhalos can be detected by the observations 
of Parkes Pulsar Timing Array (PPTA). 



Subject headings: cosmology: theory 
— dark matter 



-methods: numerical — Galaxy: structure 



Introduction 



In our Universe, dark matter halos evolve in the hierarchical fashion. Smallest microha- 
los form first, and they merge with each other to form larger halos. The size of the smallest 
microhalos is determined by the s cale of coUisional damping and free strea r aing of dark mat- 
ter particles. Analytical studies (jZybin et al.l Il999l : iHofmann et al.l l200ll : iBerezinsky et al. 
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20031 : IGreen et all 120041 : iLoeb fc Zaldarriagal l2005l : iBerezinskv et al.l 120081 ) predicted their 
mass to be 3.5 x 10~^ ~ 8.4 x 10~^Mq, though uncertainty in theory of supersymmetry 
enlarges the interval further. 



Early studies ( iBerezinskv et al.ll2003l : iDiemand et al.ll2005l : iBerezinsky et al.ll2008[ ) sug- 
gested that a significant fraction of microhalos born in early universe have survived up to 
present time, and they might be observed as the dominant sources of the a nnihilation signal. 
Thes e microhalos could enhance the annihilation signal by a factor of 2 to 5 (IBerezinsky et al. 
20031 ). who se signature rnight h ave alre ady been observed as electron a nd positron exces s 
[PAMELA (lAdriani et al.lboooh . ATIC Jchang et al.l 120081 ) . PPB-BETS JXorii et al.ll2008l ). 
and Fermi jAbdo et al.lboOQbl ) ]. 



Diemand et al.l (|2005[ ) simulated the formation of microhalos using A^-body simulations. 



They argued that the density profiles were well fitted by a single power law, p{r) oc r~'^ , with 
slope 7 in the range of 1.5 to 2 down to the radius of ~ 10~^ pc and that most of microhalos 
will survive against galactic tidal field and encounters with stars. 



On the other hand, ISpringel et al.l (120081 ) argued that the fraction of the mass in sub- 
halos in the solar neighborhood was significantly lower than that averaged over the entire 
halo. They assumed that the spatial distribution of the microhalos is the same as that of 
subhalos with mass IO^Mq, and concluded that the microhalos have a negligible impact on 
detectability. However, it is not clear whether their result for subhalos with a mass 10^M( 
can be used to estimate the distribution of halos of 10~^Mt7,. 







Since earth-mass subhalos contain no substructures by definition, their central struc- 
tures can be completely different from that of more massive halos which contain many 
substructures. However, there is no simulation of earth-mass halos with sufficient resolution 



to study the their central structure, so far. iDiemand et al.l (l2005l ) used mass resolution of 
1.2 X lO^^'^M©, which was too low to determine the central structures of microhalos. 

In this letter, we report the result of cosmological A^-body simulations with 100 times 
better mass res olution and 20 times better spatial resolution compared to that used in the 
previous work ( IDiemand et al.ll2005l ). 



2. Initial Conditions and Numerical Method 

We performed two high resolution cosmological A^-body simulations. The number of 
particles is 1024^. The size of the simulation box is 30 comoving parsecs with periodic 
boundary condition. The mass of particles is 9.43 x lO^^^M©. For the time integration, we 
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used the GreeM code (llshiyama et al.ll2009al Jbl). We used a leapfrog integrator with shared 
and adaptive time steps. The step size was determined as m\n{^Je/\ai\,e/\vi\) (minimum of 
these two values for all particles). We simulated them from z = 500 to 31. The (plummer) 
softening length e was constant in the comoving coordinate from z = 500 to z = 100, and 
constant in the physical coordinate (5 x 10~^ pc) from z = 100 to z = 31. To gen erate 



the initial particle distribution, we used the MPGRAFIC package (jPrunet et al.ll2008l ). We 
considered two models with different initial matter power spectra. 

The matter power spectrum in Mod el A includes the sh arp cutoff corresponding to dark 
matter particle with a mass of lOOGeV (iGreen et al.ll2004l ). The power spectrum of Model 
B is without cutoff. The cosn aological parameter s adopted are based on the concordance 
cosmological model [WMAPl JSpergel et al.ll2003h . Qq = 0.268, Aq = 0. 732, hn = 0.71 as = 
0.9]. However, we us ed ag = 0.8 to be close to recent observational ones (ISpergel et al.ll2007l : 



Komatsu et al.ll2009l ). We did not put any thermal v elocities in initial setup of Model A to 
avoid unphysical density fluctuations in small scales (jColm et al.ll2008l ). 



Figure [T] shows the snapshots of our simulated Universe at 2 = 31. We can see that 
there are no substructures in one halo of model A, except for caustics generated by non- 
linear growth of long- wavelength fluct uations. This structure is quite similar to what we see 



in warm dark matter simulation (e.g. iBode et al.ll200ll : iGao &: Theuna 12007), and different 



2008 
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from that of galaxy-sized or larger dark-matter halos (iSpringel et al 
2009bl : IStadel et al.l 120091 ). The difference comes simply from the initial condition 
microhalos do not contain smaller fluctuations inside. 



Our 



Results 



3.1. Structures of Microhalos 



Figure [2] shows the spherically-averaged structures of three most massive halos. We also 
performed a low resolution simulation of model A for convergence check (the number of par- 
ticles was 512'^ and the softening length was 1.0 x 10~^ pc). Two results are identical outside 
the softening radius. The central density slope for model A is considerably steeper than that 
for model B. The result of model B is similar to t hat of large-scale c o smological simulations 



(Navarro et al. 



Stadel et al 



1996 



2009 



Fukushige fc Makinol 119971 : iMoore et al.l 1 19991 : iFukushiee et all 12004 



Navarro et al.ll2010l ). The result of model A is consistent with being a 



single power law of p oc r for quite a wide range in the radius. Our result is reliable 
down to 10~^ pc, where the density reaches more than 5OOM0/pc^. 



The left-lower panel of flgure[2]shows the rough estimate of the phase-space density given 
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Fig. 1. — Top and bottom panels show the distribution of dark matter at 2; = 31 for our 
standard model (A) and model without free-streaming cutoff (B). The width of the images 
correspond to 12 comoving pc. Images in the squares are enlargements of single halos. The 
size of the squares is 0.6 comoving pc (4900AUs in physical units). 
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radius [pc] radius [pc] 

Fig. 2. — Radial profiles of three different microhalos of model A (black) and B (blue) 
at z = 31. Four panels show the density, velocity dispersion, phase-space density, and 
accumulated mass. Circles show profiles of a low resolution simulation (A512: 512^ particles, 
1.0 X 10~^ pc softening length). Two out of three profiles (middle and bottom) are vertically 
shifted downward by 1 and 2 dex. In the panel of phase-space density profiles, the two 
profiles (middle and bottom) are shifted by 2 and 4 dex. 
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Fig. 3. — Tlie differential probability distribution of the microhalo formation redshift pre- 
dicted by Press-Schechter theory. The longest wavelength correspond to lOOMpc (black), 
3k pc (green) and 30pc (r ed). These value corresponds to the box length of the simulations 
by lDiemand et al.l (|2005[ ) (3kpc) and this work (30pc). 



by where (f^) is the velocity dispersion. The slope is close to —2.25 for r < 10^'^ 

pc. Because of Liouville's theorem, this phase space density cannot exceed the initial value 
at the time of decoupling, which is ~ lO^^M0pc~^(km/s)~'^. Thus, the gravitational collapse 
of the cusp stops at the radius at which the phase space density reaches this limit. This 
radius and density there are ~ 10~^pc and ~ 2 x lO'*M0pc~^. We, therefore, can safely 
conclude the density profile of microhalos is 



p(r) = Pc{r/rc) 



-1.5 



for 10 ^pc > r > Tc 



(1) 



and p(r) ~ pc for r < r^. This profile is co nsistent with the r esults of early studies of the 
collapse of self-gravitating gaseous spheres (ISuto fc Silkl Il988l ) and rec ent high resolutio n 
simulations of cold collapse ( Nipoti et al. 20061 ) and warm dark matter ( Colin et al. 2008 ). 



The average density of these earth- mass halos is affected by density fluctuations of larger 
scales. However, our simulations do not c ontain large scale fluctuations. As seen in Figure [31 
Press-Schechter (jPress fc Schechterlll974j ) theory predicts that the average formation epoch 
of microhalos is earlier if we take into account larger scale fluctuations. Since the average 
density of a microhalo should be proportional to the cube of its formation redshift, most 
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microhalos should have higher average density than these in our calculation. For simplicity, 
we use the profile of equation ([1]) to estimate the chance of survival of these halos. 



3.2. Tidal Cutoff and Encounters with Stars 



The tidal cutoff radius of microhalo with the density profile of equation ([T]) due to 
galactic tide is expressed as rmax = 0.082(i?/10kpc)^/^pc, where R is the distance from the 
galactic center, if we assume an isothermal halo with the rotation velocity of 220 km/s. We 
can see that these halos can easily survive at the distance Ikpc, and even at the distance 
O.lkpc, the central 10~^ pc of the halo would still survive. 



Can encounters with stars destroy the internal cusps of microhalos (IBerezinskv et al. 



2003 



2007 



2006 



Angus fc Zhao 120071 : ICoerdt et all 120071 IZhao et all 120071 ICreen fc Goodwin 



Schneider et al.l 120101 )? Since the surface number density of disk stars is about 100 



pc ^ at the solar neighborhood, the closest distance a halo with the orbital period of the 
order of 100 Myrs can approach to a star in the Hubble time is 10~^ pc. 

FigureHlshows the result of encounter with a IMq star moving at 200kms~^. We selected 
the most massive microhalo at ^ = 31 from the simulation of model A (the cutoff radius 
is two times virial value), and added velocities on each particle according to the impulsive 
approximation. Then, we simulated the evolution of the microhalo for 27Myr after the 
encounter. We can see that the regions of radius ~ 10~^ pc survives after encounters with 
the impact parameters of 0.02 pes. Theoretically, the cutoff radius due to encounter for 
the power-law cusp with slope —1.5 is given by r oc b^^^^, where b is the impact parameter. 
Therefore, the complete disruption requires an encounter with 6 = 5 x lO^^pc. Such close 
encounters are expected only at the very central regions of the galaxy (20pc or less from the 
center) . 



If we extrapolate the subhalo mass function dn/dm oc (IBerezinsky et al.ll2003l ) to 
the mass of microhalos, the number of microhalos in our Galaxy is ~ 5 x 10^^ (we assumed 
the mass function dn/dm ~ 5 x 10^°m~^). The number density of the microhalos near the 
solar neighborhood is ~ 500pc~^, if the distribution of the microhalos is the same as that of 
the background dark matter. This extrapolation is justified by the fact that the microhalos 
survive under tides from our Milky Way Galaxy. Note that the above estimate takes into 
account the effect of merging between small halos and accretion to somewhat larger halo 
since we extrapolated the mass function obtained by iV-body simulation. 



Fig. 4. — Snapshots and radial density profiles after encounters with stars. Panel a (top- 
left) shows the snapshot of a halo in model A at z = 31. Only particles within the radius 
of 0.037 pc from the center of the halo are shown. Panel b (top-middle), c (top-right), d 
(bottom-left), and e (bottom- middle) show the distribution of particles at 27 Myrs after 
the encounter with a solar-mass star with the impact parameters 5 x 10~^pc, 1 x 10~^pc, 
5 X 10~^pc, and 1 x 10~^pc. The black regions in the center are caused by the initial outer 
cutoff of the particle distribution and are not real. The upper right box in each panel shows 
the central region of the microhalo. The central regions are located at the center (Panel 
a and b), and the positions indicated by arrows (Panel c and d). In panel e, the central 
region is destroyed completely. Panel f (bottom-right) shows the density profiles after the 
encounter. 
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3.3. The Gamma-ray Luminosity of a Microhalo 

The gamma-ray flux of a microhalo is given by 

N lm)\ 

d \n{r^^^/r^)pj^, (3) 



2ml 



2 



30 3 X 10-26cm3s-i VlOOGeV/ V / 

' ' 1 / 4 g]^ _|_ Jli^^ — J photons ■ cm~^ ■ s~^, 



2 X lO^Mopc-V VlO^^Pc/ v' rc/10-5pc 

where A^^ is the number of emitted photon per annihilation, is the mass of dark matter 
particle, {av) is the interaction cross section of dark matter, and r^ax is the outer cutoff 
radius of a microhalo at which the slope of the density becomes steeper than —1.5. 

Mass loss due to the galactic tide or encounters with stars changes rmax of microhalos, 
but does not affect pc or '^c except in the cases of extremely close encounters, as we can see 
in flgure 3. Therefore, they have relatively minor effect on the gamma-ray luminosity. After 
99% of the mass is lost from a microhalo, it still retains nearly 50% of the luminosity. Thus, 
we can assume that the tidal stripping has practically no effect on the gamma-ray luminosity 
of our earth-mass halos . These results supports the predictions of early analytical studies 



flBerezinskv et all 120081 1. 



Here, we consider how the gamma-ray flux of the microhalo depends on its formation 
epoch. The average density of a halo pave reflects the cosmic density at its formation time 



( Bullock et al. 2001 ). Since the cosmic density is proportional to (1 + 2f)^, pavc oc ^y+f^) ' 



where and zq is the formation redshift and the typical formation redshift of microhalos in 

— 1/3 

our simulations. From the conservation of the mass, we can derive rmax pave • We can 
rewrite p^r^ in equation (jll) as Pave'^max using equation ([T]). Therefore, 

Here, we drop the dependency of the formation redshift in the logarithmic term in equation 
since it is rather weak (vr oc rmax'^). 



The average boost factor due to the formation epoch is determined as 

-dzf, (6) 



dP{zi) F^{zi) 



dzi F^(40) 
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where P{z{) is the distribution function of the formation epoch. For the PS function in Figure 
[3l its value is ~ 1.6. We apply this boost factor to estimate the gamma-ray luminosity. 



4. Discussions and Summary 
4.1. Gamma-ray Signal from Microhalos 



There are many works o n whether or not subhalos and microhalos can be observed via 



annihilation gamma-rav fe.g.lOda et al 



2005 



Koushiappasll2006uAndo et al.ll2008l : iLee et al. 



20091 : iKamionkowski et al.l l2010l : ISchneider et al.l |2010[ ) . We made the all sky map of the 
gamma-ray annihilation signal in our Galaxy based on the new profile of the microhalo. The 
observer locates at 8.5 kpc from the center of the halo along its long axis. We consider the 
emissions from microhalos only. We assume the spatial distribution of microhalos follows 
the mass distribution of a galaxy sized LCDM halo selected from our previous simulation 
(jlshiyama et al.ll2009al ). Thus, our all-sky map naturally include all subhalos with the mass 
larger than numerical limit (~ 10^ Mq). In addition, destruction due to encounters with 
stars are taken into account. We assumed a exponential disk with the surface density, 
S(r) ~ 1000 exp(— r/3kpc)pc~^, where r is the distance from the galactic center. From the 
position of an ensemble of microhalos, the impact parameter b of encounters with stars are 
calculated. Then, we use the cutoff radius r^ut = 0.0176^/^^pc and reduce the gamma-ray 
luminosity. 

In order to include the luminosities of nearby microhalos, we placed randomly 2000 
microhalos less than Ipc from the observer. The distribution of formation redshift and the 
gamma-ray flux of microhalos are given by Figure [3] and equation ([5]). 

Figure O shows the all sky map. Many nearby microhalos are observed as pointlike 
sources, since their angular size is around 1'. Therefore, the nearby microhalos are promising 
sources of gamma-ray signal of the annihilation of dark matter particles. Their distance is 
around 10^ AU and velocity is of the order of 200km/s. The proper motion is as large as 
0.2 deg/year. They might be too dim to be observed by Fermi. Sommerfeld effect could 
boost the interaction c ross section of dark matter, and enhance by orders of magnitude the 
gamma-ray luminosity f Kuhlen et al. 20091 ). These microhalos might be good targets of next 
generation Cherenkov telescope, such as ctaEI. 



Bright microhalos enhance luminosities of all subhalos (including known Milky Way 



^http: / / www.cta-observatory.org/ 



-1 



log Flux [photons m ^ sr ^ ] 




Fig. 5. — All sky map of the gamma-ray annihilation signal of our Galaxy. Pointlike 
sources indicated by arrows show the contributions from microhalos less than Ipc from the 
observer and with flux larger than 0.03 photons ■ m^^ ■ s^^ • sr~^. Here, we assume that the 
number of photons per annihilation is Nj = 30, the interaction cross section of dark matter 



is {av) = 3.0 X 10 



-26 



■H —1 

cm s , 



and the mass of dark matter particle is rriy = lOOGeV. The flux 



of the integrated emission of microhalos in the galactic center is 1.40 photons ■ m~^ ■ s~^ ■ sr~^. 
The flux of the most brightest microhalo is 0.069 photons ■ m~^ ■ s~^ ■ sr~^. Flux needed for a 
lOo" detection for Fermi is ~ 0.4 x 10~®photons ■ cm~^ ■ s~ ^ for low Galactic latitude and is 



~ 0.1 X 10 ^photons ■ cm ^ ■ s ^ for high Galactic latitude ( lAbdo et al.ll2009al ). The angular 
resolution of Fermi is about 0.1 degree for energies above IGeV. 
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satellite) since subhalos contain a number of microhalos. The brightest gamma-ray source 
is the integrated emission of microhalos in th e galactic center, but it is much more extended 
compared to the result of previous studies (ISpringel et al.l l2008l ) . These sources may be 
also good targets for observation. However, baryons would have significant effect on the 
distribution of dark matter at the center of the galaxy. If baryons dis rupt the central cusp of 
the main halo ( iMashchenko et al.ll2006l . l2008l : iGovernato et al.ll2010l ). gamma-ray flux from 
the galactic center might be muc h smaller. Observa tional data suggests that dwarf galaxies 
have constant density cores (e.g. iGentile et al.ll2004j ). 



4.2. Perturbations on Millisecond Pulsars 



Pulsar timing measurements might be used to detect microhalos. Since many of known 
millisecond pulsars (MSPs) are in the direction of the galactic center, the number density of 
microhalos around MSPs is higher than that of solar neighbourhood. 

Here, we estimate the residual of the time of arrival of a MSP d ue to the gravitational 
attraction from its nearest neighbor, following ISeto fc Coorayl ( 120071 ). The residual due to 
the constant acceleration a in the period of t is given by 



AT 



2^ 



(7) 



where c is the speed of light, and a is given by 



GM 

« = -^, (8) 

where G, M and R are the gravitational constant, the mass of the microhalo, and the distance 
between the microhalo and the MSP. So we have 



AT = 10.8 



R 



2xl04AUy VIO-6M0 



M 



t 



20yr 



ns. 



(9) 



If we assume that a microhalo lost 90% of mass by encounters with stars, the residual is 
about ~ 2ns in twenty years. The change of the distance in ten years is around 400AU. 
Since the change in the acceleration is proportional to the third power of the distance, one in 
ten MSPs would show the change in the accel eration ten times l arger than that of a typical 
MSP, which should be detectable with PPTA ( IManchesterl l2008l ) after modest improvement 
in the timing accuracy. 

Note t hat the mass loss of two weak encounters can be larger than that of a strong 
encounter (lAngus fc Zhaol 120071 ) . For simplicity, we used the mass loss due to a strong 
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encounter as estimated in subsection 13.21 However, the mass loss of microhalos in the 
opossite direction to the galactic center is less than that near Sun. MSPs in the direction 
should show larger timing residual. 



4.3. Summary 

We found that earth-mass microhalos have steep central density cusps. Their central 
regions is not disrupted by tidal forces from the parent galaxies or stars and they survive to 
the present time. Our result is different from the re cent claims that smal l-scale structure have 



a negligible impact on dark matter detectability ( ISpringel et al.ll2008l ). They considerably 
underestimate annihilation signals because they assumed the density profile for the smallest 
microhalos was shallower than that we found with high resolution simulations. 

We are grateful to Veniamin Berezinsky for helpful discussions. Numerical computations 
were carried out on Cray XT4 at Center for Computational Astrophysics, CfCA, of National 
Astronomical Observatory of Japan. T.I. is financially supported by Research Fellowship of 
the Japan Society for the Promotion of Science (JSPS) for Young Scientists. This research is 
partially supported by the Special Coordination Fund for Promoting Science and Technology 
(GRAPE-DR project). Ministry of Education, Culture, Sports, Science and Technology, 
Japan. 
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